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We demonstrate how a time- dependent analysis of the decays D° — > 
K^K^it^ can be used to determine the D° — D° mixing parameter y with 
a precision that is competitive with established methods. The proposed 
analysis is an inclusive study which makes use of the measurements of 
^ | the coherence factor and mean strong phase difference for these decays 

recently performed by CLEO. 
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1 Introduction 



In the last few years measurements have been performed at e + e~ colliders and the 
Tevatron which, when taken together, reveal the presence of mixing in the D° — 
D° system. These measurements have relied on a variety of analysis strategies, 
including the time-dependent study of D° — > K + n~ [H 12], D° — > K + n^n° and 
D° — > K + 7t~7t + 7t~ [3] decays, the measurement of the D° lifetime in decays to 
CP-eigenstates [4L and the amplitude analysis of self-conjugate final states, such 
as Kg7i + 7i~ [SJ E] □• The mixing is parameterised through the dimensionless variables 
x and y, which determine the evolution of a D° into a D° meson through off-shell and 
on-shell intermediate states respectively The current values for these parameters, 
averaged over all analyses, are x = (0.63loio)% an d V = (0.75 ± 0.12)% [7]. More 
sensitive individual measurements of x and y are needed, both to test the consistency 
of the overall picture, and to search for evidence of CP-violation in D° — D° mixing, 
where the existing constraints are very weak [?]. As CP- violation in the D° — D° 
system is generally considered to appear at a negligibly small level in the Standard 
Model [8j, any observation of such an effect would constitute a clear signal of New 
Physics. 

In this Letter we propose to study D° — D° mixing through a time-dependent 
inclusive analysis of the decays D° — > KgK~7t + and D° — > KgK + 7r~, in a manner 
similar to that of a 'wrong sign' D° — > K + n~ analysis. Whereas the Cabibbo favoured 
(CF) and doubly-Cabibbo suppressed (DCS) amplitudes in the Kti final state give 
rise to time integrated branching ratios that differ by a factor of ~ 300, for the 
singly- Cabibbo suppressed (SCS) modes there is approximate equality between the 
two final states, with the branching ratios of D° — > KgK~7i~ and D° — > KgK + 7i + 
being (3.5 ± 0.5) x 10~ 3 and (2.6 ± 0.5) x 10~ 3 respectively [9]. This attribute gives 
rise to several interesting and distinctive features. 

Section |2] provides details on the method, and the sensitivity of the analysis is 
assessed in Section [3J Conclusions are presented in Section HI Throughout the dis- 
cussion, and by way of analogy with the CF/DCS strategy, the D° decay to the K + 
(K~) final state is designated as 'wrong sign' ('right sign'), and it is assumed that 
the direct decay of the D° to the wrong sign (right sign) final state proceeds by the 
suppressed (favoured) amplitude. 



2 Measuring y with the decays D — > K s K T n 

Adopting the notation in Ref. [10], and assuming no CP violation, the D mass 
eigenstates, | Di^), can be written in terms of the flavour eigenstates | D°) and | D°) 

*Unless stated otherwise, the charge-conjugate mode is implicit throughout this Letter. 



1 



as follows: 



\D 1 ) = 


\D°) + 


\D°), 


\D 2 ) = 


D ) — 


\D°). 



The time evolution of the physical states is given by 

| Djit)} = e - iM ^ T ^ | Dj(t = 0)), (2) 

where Mj and Tj are the masses and widths of the mass eigenstates. 
We are interested in the time evolution of a sample of initially pure 

I D } mesons 

decaying to the wrong sign final state | /). The decay can occur directly from the 
D° or via mixing followed by decay from the D°. The total decay amplitude is 
the sum of the amplitude from the direct decay (the suppressed amplitude) and the 
amplitude from | D°) to the final state (the favoured amplitude). The final state 
under consideration is multi-body (n > 3) and hence the amplitudes are functions 
of the final state kinematics, p. We write the suppressed amplitude as A(p) = 
(/ p | H | D°) and define A 2 = J ^.(p)^4(p)dp. The favoured amplitude is written 
as B(p) = (f p | H | D°), with B 2 = J B(p)B(p)dp. The analysis is inclusive in the 
sense that there is no attempt to distinguish the intermediate resonances contributing 
to | /). Using the definition of the 'coherence factor' Rd [TT] we can write that 
ABR D e- iS ° = J A(p)B(p)dp, where S D is the average strong phase difference between 
the suppressed and favoured amplitudes. 

Defining the following variables in the usual way: 

(3) 
(4) 
(5) 
(6) 
(7) 



M = 


^(M 1 + M 2 ), 


AM = 


Mi -U.. 


r = 


^(Ti + r a ), 


Ar = 


r 2 -r 1; 




AM 


X = 


r ' 




Ar 


y = 


2T' 


we can write the time evolution of the state D°(t)) which is a pure sample of 


at t — 0, using expressions (pQ) and ^ 


as 


1 D°(t)} = f + (t) 


\D°) + U(t)\D°), 



(9) 
where 

f + (t) = e - iMt -^cos(^(x-iy)), (10) 
=e- iMt -? r Hsm(% (x + iyj). (11) 
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Hence the total decay amplitude to the wrong sign final state | / p ) with kinematics 
p is then 

(/ p | H | D\t)) = A(p)f + (t) + B(p)f-(t), (12) 
and the total decay rate to the wrong sign state is then given by 

r[D°^f] = j[A(p)f + (t)+B(p)ut)}x[^f4tr+w^f-(tr}d P , 



(13) 



= e 



-rt 



A 



B 2 



^-(e yrt + e- yTt + 2 cos(xFt)) + ~^{e yVt + e~ yVt - 2 cos(aTt)) 



(e~ i6F [e yTt - e~ yTt - 2i sm(xTt)] + e i5p [e yFt - e~ yVt + 2i sm(xFt)] ) 

Let us define td = A/B and y' = y cos Sd — xsm5o and specialise to the case 
D° — > KgK^n^. Under the assumption that i< 1 and y <C 1, the decay rate to the 
wrong sign state up to orders of x 2 and y 2 becomes 



T[D° ->■ K° s K + tt- 



-rt 



+ 



K° s Ktt K°Kn , 



-(rtf 



(14) 



where labels indicate the KgK T ir ± specific quantities. We may also consider the rate 
to the right sign state, | f p ). In this case, the decay amplitude is given by 



</ p | H | D°(t)) = A(p)f + (t) + B(p)f-(t). 



(15) 



Once more specialising to the case D° — > KgK T n ± and assuming x C 1 and y <C 1 
we obtain 



+ 



.^S K *\2 „2 



l-(r^") 2 U 2 



■(Tt? 



• (16) 



It is instructive to compare with the D° — > K^-k^ case. There « 1 in con- 
trast to the SCS decays where r^ sKw ~ 1. Furthermore, being a two body decay the 
coherence factor reduces to unity, and is the strong phase difference between the 
suppressed and favoured amplitudes. The familiar wrong-sign Kn mixing expression 
is then obtained: 



T[D° -> K + tt-\ = 



(17) 
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whereas the right-sign rate has negligible deviation from an exponential decay. 

Several observations can be made when assessing the potential of D° — > K^K^tx^ 1 
in a mixing analysis: 

• Recently the CLEO collaboration has presented the first measurements of the 
coherence factor and mean strong phase difference in D° — > K ( $K ± 7[ T de- 

cays pj]. Preliminary results of Rp = 0.73 ±0.09 and 5p = (8.2 ± 15.2)° 
are reported. The relatively large value of the coherence factor means that the 
interference term in y' K o KlT in expressions (JT4l) and (Tl6|) receives little dilution 
from the presence of the intermediate resonances. The measured value can be 
used as a constraint in a mixing analysis. Furthermore, the measured value of 
fiK s Kn k e usec [ ^ re i a te y' K o K7T and y. 

• The fact that r^, sKlT » has several consequences when comparing D° — > 
K'gK Tr K ± with D° — > K T n ± as a mixing channel. Firstly in the D° — > KgK T n ± 
case both the wrong sign and right sign final states have sensitivity. On the 
other hand, the relative size of the term in Tte~ rt is small, giving the mode 
lower event-by-event sensitivity to y' K o K7T than the Kir analysis has to y' K7T ; 

furthermore the (rt) 2 e~ r * term is negligible. Therefore D° — > KgK^Tt is 
effectively sensitive to y' K o Kn alone, provided the sample is sufficiently large to 
fit the deviation from the purely exponential decay. 

• Experimentally, it is easier to obtain a high purity sample when selecting D° — > 
KgK 11 ^ than D° — > K + n~ . This is because of the higher branching ratios, 
and because D° — > K + n~ receives significant contamination from D° — > K + n~ 
decays in events where a random pion of the wrong charge is mistaken as the 
'slow pion' when reconstructing the D* + — > D tt + decay chain. 

3 Sensitivity studies 

The sensitivity of a mixing measurement using D° — > K^K^it^ decays is assessed 
via a simple Monte Carlo study. Comparisons are made between a published wrong 
sign D° — y K + tt~ analysis performed at BABAR [2] with 384 fb _1 and a possible 
measurement using D° — > K^K^-k^ exploiting the same dataset. 

The true decay time distributions of simulated signal data are generated according 
to the distributions given in expressions f|T4]) . ( ITS]) and f lT7|) . where the values of the 
constants required for the simulation are given in Table [1] and are taken from [7] , [9] 
and [LJ. 

The detector resolution is simulated by assigning each event a decay time un- 
certainty according to a Landau distribution with most probable value 0.16 ps and 
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Table 1: Input values used in the simulation study. 



Parameter Value 
t(£> u ) 0.410 ps 

y 0.0075 
x 0.0063 
0.057 
6%* 22° 



rn s 0.88 



tit** go 

i?^^ 0.73 



<7 = 0.04 ps. Values of the decay time uncertainty are restricted to those less than 
0.5 ps. This description of the decay time uncertainty is similar to that described in 
[2J. The generated decay time is then smeared by drawing a random value from a 
Gaussian with a width of the event decay time uncertainty. 

Three categories of background events are considered. We assume that 75% of 
the background events are true D° decays where a random slow pion has been used 
to tag the event incorrectly. A further 15% of the background events are assumed to 
be combinatoric and have a true decay time of 0. The remainder of the background 
is assumed to be from partially reconstructed charm decays and the time distribution 
is generated according to an exponential decay time with lifetime = 0.2 ps. The 
decay time uncertainty is assigned using the same distribution as for signal events. 
The resulting total decay time distribution for background events is similar to that 
presented in [2J. 

In 468.5 ftr 1 of data accumulated at BABAR, 540,000 tagged signal D° -»■ 
KgiT + 7c~ events are observed with 98.5 % purity [6]. Assuming the branching frac- 
tions ratios in [9] we extrapolate that the same dataset used for the D° — » K^-k^ 
mixing measurement (384 fb _1 ) should yield 80,000 D° — > K^K^n^ signal events 
with purity better than 98%. The signal events are split between the wrong and right 
sign according to their measured branching fractions. To compare the sensitivity 
with existing measurements using D° — > K + tt~ decays we generate a similar Monte 
Carlo assuming 4000 signal events with a purity of 50% as found in [SJ. The decay 
time distributions of the simulated wrong sign samples are shown in Figure [1] where 
the components due to direct decay, interference, mixing, and background are shown 
separately. The difference in contribution of the interference term and background 
for the two final states is clearly visible. 

To determine the sensitivity to y' K o K , 1000 independent simulated samples of 
events are generated and fit. A maximum likelihood fit is performed on each sample 
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where the only floating parameter is y' K o K and all other constants are fixed to the 
values given in Table [1] and those describing background. The likelihood ignores 
the second order terms in x 2 and y 2 , as there are insufficient data to extract these 
parameters also. The RMS spread on the returned value of y' K o K , corresponding to 
the expected statistical sensitivity, is 0.0054. However the mean of the fitted values 
is 0.0014 lower than the input, indicating a ~20% bias. The source of this bias has 
been studied. It is not a consequence of ignoring the second order terms; rather it is 
found to be due to an intrinsic limitation in correctly extracting the fraction of events 
distributed as te _i//T in the presence of events distributed as e~ t//T when the fraction 

is very small. The magnitude of the bias depends on the relative size of r^ sKn and 
y' K o Kn , but is found to be approximately constant within the experimentally allowed 
ranges of these parameters. It is assumed that this effect can be corrected for in the 
measurement. 

There is a source of systematic error due to the current 12% uncertainty in the 
knowledge of the coherence factor. This is found to induce an uncertainty of 0.0010 
in the fitted value of y' K o K ■ Finally, it should be noted that as the time distribution 

of both wrong and right sign D° decays is sensitive to y' K o Kn , the mistag background 
will also contain events from the interference of mixing and the suppressed decay. 
As the purity of the sample is high, approximating the mistag decay time by an 
exponential causes negligible bias. The total uncertainty on a y' K o Kn measurement 
using 80,000 signal decays is therefore predicted to be 0.0055 taking into account 
statistical uncertainties and the present knowledge of the coherence factor. 

The measured value of y' K o K can be related to the mixing parameter y with 

K° Ktt 

knowledge of x and the mean strong phase difference, S D S . The sensitivity to y 
using the known values and uncertainties of these parameters is 0.0059. 



6 



A fit study is also carried out on the generated D° — > K + tt~ events. In this case 
the second order term in the lifetime distribution is not ignored and both y' K% and 
x' 2 are floating parameters. The correlation between the fit parameters is found to 
be -0.94 and the sensitivity to y' K7T and x' 2 is 0.0046 and 0.00026 respectively. These 
values are very similar to those reported in [2], validating our procedure. We can 
translate the result for y' Kn into a measurement of y and find a total uncertainty of 
0.0052. 

4 Conclusions 

A time dependent inclusive analysis of D° — > K^K^n^ decays appears a promising 
method to improve our knowledge of the D° — D° mixing parameter y. Although the 
relative size of the interference effects is significantly smaller in these modes than in 
D° — > K + 7i~, compensation is provided by the higher branching ratio and expected 
high sample purity. Recent preliminary results from CLEO indicate that the dilution 
arising from intermediate resonances to the interference term is not severe. The 
analysis can be generalised to look for evidence for CP-violating effects in D° — D° 
mixing. 

A sensitivity study based on existing 5-factory publications shows that the preci- 
sion obtainable on y with D° — > K^K^Tt^ decays is similar to that of the D° — > K + n~ 
analysis. Experimentally, careful checks of the proper time determination will be re- 
quired to ensure that systematic effects do not bias the measurement. Validation 
could come from demonstrating a good understanding of the proper time distribution 
of the mixing-free decays of non-.D° mesons to final states of similar topology. 
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